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Abstract 


Accurate  ground  moving  target  indication  (GMTI)  and  target  parameter  estimation 
can  be  achieved  only  after  sufficient  suppression  of  interfering  stationary  clutter,  par¬ 
ticularly  for  space-based  SARs  with  typically  small  exo-clutter  regions.  In  its  simplest 
form,  this  is  accomplished  using  two  radar  receiver  channels,  such  as  the  dual  receive 
antenna  mode  of  RADARSAT-2’s  Moving  Object  Detection  Experiment  (MODEX). 
In  this  mode  of  operation,  the  full  antenna  is  broken  up  into  two  sub-apertures  with 
two  parallel  receivers  to  create  two  independent  phase  centres.  It  is  well  known,  how¬ 
ever,  that  a  two-aperture  approach  to  GMTI  is  sub-optimum  and  that  target  param¬ 
eter  estimation  is  often  compromised  by  clutter  interference  or  poor  signal-to-clutter 
ratios.  Two  degrees-of-freedom  are  simply  not  enough  to  simultaneously  suppress 
the  clutter  and  to  accurately  estimate  the  target’s  properties,  such  as  velocity  and 
location. 

The  investigation,  described  in  this  Technical  Memorandum,  explores  several  con¬ 
cepts  of  increasing  the  spatial  diversity  for  RADARSAT-2,  which  allows  the  two- 
channel  SAR  system  to  operate  like  a  three  or  four  channel  radar.  Owing  to  the 
very  flexible  programming  capabilities  of  the  RADARSAT-2  antenna,  this  can  either 
be  accomplished  by  the  toggling  of  the  transmitter  between  subsequent  pulses  or  via 
clever  transmitter /receiver  excitation  schemes.  A  trade-off  analysis  between  number 
of  channels,  phase  centre  separations,  and  PRF  limitations  is  presented  for  a  system 
based  on  RADARSAT-2  MODEX  parameters. 


Resume 


La  precision  des  indications  de  cible  terrestre  mobile  (GMTI)  et  des  estimations  de 
parametres  des  cibles  depend  uniquement  d’une  suppression  suffisante  du  clutter  fixe, 
surtout  pour  les  ROS  spatiaux  a  regions  d’exoclutter  typiquement  reduites.  Sous  sa 
forme  la  plus  simple,  ce  resultat  s’obtient  au  moyen  de  deux  canaux  de  recepteur  ra¬ 
dar,  par  exemple  dans  lc  mode  a  antenne  de  reception  double  utilise  pour  l’experience 
de  detection  des  objets  mobiles  (MODEX)  de  RADARSAT-2.  Dans  ce  mode  de  fonc- 
tionnement,  l’antenne  complete  est  subdivisee  en  deux  sous-ouvertures  associees  a 
deux  recepteurs  paralleles  creant  deux  centres  de  phase  independants.  On  salt  bien, 
toutefois,  que  la  technique  des  deux  ouvertures  n’est  pas  optimale  pour  les  applica¬ 
tions  GMTI  et  quo  l’estimation  des  parametres  des  cibles  est  souvent  compromise  par 
le  brouillage  du  au  clutter  ou  par  de  mauvais  rapports  signal/clutter.  Deux  degres 
de  liberte  ne  suffisent  tout  simplcment  pas  pour,  simultanement,  supprimer  le  clutter 
et  estimer  avec  precision  les  proprietes  des  cibles,  comme  la  vitesse  et  la  position. 
L ’etude  decrite  dans  lc  present  document  technique  examine  plusieurs  aspects  de 
l’accroissement  de  la  diversite  spatiale  de  RADARSAT-2,  qui  perrnet  de  faire  fonc- 
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tionner  le  systeme  SAR  a  deux  canaux  comme  un  radar  a  trois  ou  a  quatre  canaux. 
En  raison  dc  la  tres  grande  souplesse  de  programmation  de  l’antenne  RADARSAT-2, 
ce  resultat  pent  s’obtenir  soit  par  le  basculement  de  Pemettenr  entre  des  impulsions 
subsequentes,  soit  par  d’ingenieux  systemes  d’excitation  d’emetteur  et  de  receptenr. 
Une  analyse  de  compromis  entre  le  nombre  de  canaux,  l’espacement  des  centres  de 
phase  et  lcs  limites  de  frequence  de  repetition  des  impulsions  (FRI)  est  presentee  pour 
un  systeme  base  sur  les  parametres  MODEX  de  RADARSAT-2. 


ii 


DRDC  Ottawa  TM  2006-209 


Executive  summary 


A  Simulation  Study  of  Multi-Channel  RADARSAT-2 
GMTI 

S.  Chiu;  DRDC  Ottawa  TM  2006-209;  Defence  R&D  Canada  -  Ottawa;  November 

2006. 

Background:  Canada’s  RADARSAT-2  spaceborne  SAR  will  be  launched  in  2007. 
This  SAR  can  be  operated  in  an  experimental  ground  moving  target  indication 
(GMTI)  mode  called  MODEX  for  Moving  Object  Detection  Experiment.  In  this 
mode  of  operation,  the  radar  antenna  is  split  into  two  sub- apertures,  each  with 
its  own  receiving  channel,  which  provide  scene-motion  sensitivity  and  thus  permit 
moving  targets  to  be  detected  and  their  velocity  parameters  estimated.  However, 
a  severe  GMTI  limitation  exists  for  the  two-aperture  radar  system  clue  to  contam¬ 
ination  of  moving  target  signals  by  clutter.  This  limitation  has  motivated  the  ex¬ 
ploration  of  multi-channel  concepts  using  aperture  toggling  or  switching  techniques, 
which  allow  a  radar  with  two  physical  apertures  to  behave  like  a  multi-aperture  sys¬ 
tem  by  generating  virtual  phase  centres.  Three  multi-channel  conceptual  modes  for 
the  RADARSAT-2  type  imaging  radar  are  proposed,  examined  and  tested  in  this 
Technical  Memorandum. 

Principal  results: 

Multi-channel  concepts 

Three  multi-phase  centre  modes  for  the  RADARSAT-2  MODEX  are  proposed.  The 
first  is  a  three-channel  togglc-transmit  mode,  in  which  the  radar  alternately  transmits 
using  the  fore  and  aft  sub-apertures  and  receives  echoes  with  both  sub-apertures 
simultaneously.  The  second  and  third  multi-channel  schemes  are  three-channel  and 
four-channel  togglc-receive  modes,  respectively.  In  these  modes,  the  radar  transmits 
using  the  full  aperture  but  receives  echoes  with  different  sub-sections  of  the  aperture 
on  alternating  pulses.  The  three-channel  togglc-transmit  mode  has  the  same  along- 
track  baseline  as  the  original  two-aperture  system.  However,  for  both  of  the  three- 
channel  and  four-channel  toggle-receive  modes  the  baseline  is  half  of  that  for  the  two- 
aperture  system.  A  shorter  along-track  baseline  will  usually  lead  to  lower  sensitivity 
to  slowly  moving  targets.  On  the  other  hand,  a  higher  degree  of  freedom  (i.e.  the 
number  of  available  channels)  may  improve  the  detection  performance,  which  may 
compensate  for  this  reduced  baseline. 
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Beamwidth  and  sidelobe  level  trade-offs 


Different  combinations  of  transmit  /receive  aperture  sizes  yield  different  two-way  an¬ 
tenna  patterns  with  the  accompanying  variations  in  sidelobe  levels  and  beamwidths. 
The  widest  beamwidth  (but  not  the  highest  sidelobe  levels)  is  generated  using  half  an¬ 
tenna  transmit  and  half  antenna  receive  (i.e.  the  toggle-transmit  mode).  The  antenna 
pattern  generated  for  the  full-antenna  transmit  and  quarter  antenna  receive  (i.e.  the 
four-channel  togglc-receive  mode)  does  not  differ  significantly  from  the  full-antenna 
transmit  and  half-antenna  receive  case  (or  the  two-aperture  receive  mode). 

Improvement  in  ATI  phase  estimation  and  moving  target  detection 

Multi-channel  GMTI  performance  using  the  simple  DPCA  detection  metric  is  found 
to  improve  as  a  function  of  the  number  of  available  channels  for  a  constant  separation 
of  phase  centres.  Similarly,  the  estimation  of  the  along-track  interferometric  phase, 
which  is  proportional  to  the  radial  velocity  of  a  moving  target,  is  shown  to  improve 
in  accuracy  with  an  increasing  number  of  channels.  These  results  apply  only  to  the 
constant  baseline  case.  Increasing  the  number  of  channels  at  the  expense  of  reducing 
the  baseline  may  not  lead  to  an  overall  improvement  in  GMTI  performance. 

Multi-channel  clutter  cancellation 

The  key  question  concerning  the  generation  of  three  or  four  phase  centres  from  a  two- 
channel  SAR  system  is,  ”  Can  these  virtual  phase  centres  generated  using  aperture 
toggling  be  effectively  used  to  suppress  the  stationary  background  clutter?”  Simula¬ 
tion  experiments  which  mimic  the  three  proposed  multi-channel  modes  clearly  show 
that  the  stationary  clutter  can  be  effectively  cancelled  using  phase  centres  generated 
via  sub-aperture  toggling  in  such  a  way  that  the  interferometric  phase  of  a  moving 
target  can  be  more  accurately  estimated. 

Significance  of  results:  The  feasibility  of  applying  the  aperture  toggling  or  switch¬ 
ing  concept  to  the  RADARSAT-2  multi-channel  MODEX  is  demonstrated  here  for  the 
first  time  using  simulation.  Although  the  application  of  this  concept  to  multi-channel 
GMTI  is  not  new,  the  ” proof’  of  the  concept  is  presented  here.  This  conclusive  result 
leads  to  the  real  possibility  of  RADARSAT-2  multi-channel  MODEX  modes:  these 
were  not  considered  possible  in  the  initial  designs  of  the  MODEX. 

Future  work:  Future  work  should  include  full-scale  simulation  experiments  using 
SBRMTISIM  (the  Space-Based  Radar  Moving  Target  Indication  Simulator)  to  inves¬ 
tigate  the  advantages  and  disadvantages  of  each  of  the  proposed  modes. 
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Sommaire 


A  Simulation  Study  of  Multi-Channel  RADARSAT-2 
GMTI 

S.  Chiu;  DRDC  Ottawa  TM  2006-209;  R  &  D  pour  la  defense  Canada  -  Ottawa; 
novembre  2006. 

Contexte:  Le  SAR  spatial  RADARSAT-2  du  Canada  sera  lance  en  2007.  Ce  SAR 
peut  fonctionner  dans  un  mode  experimental  d’indication  de  ciblc  terrestre  mobile 
(GMTI)  appele  MODEX,  pour  ”  Moving  Object  Detection  Experiment  ”  (experience 
de  detection  des  objets  mobiles).  Dans  ce  mode  de  fonctionnement,  l’antenne  radar  est 
subdivisee  en  deux  sous-ouvertures,  dotees  chacune  de  leur  propre  canal  de  reception, 
qui  procurent  une  sensibilite  an  mouvement  des  scenes  et  permettent  ainsi  de  detecter 
les  cibles  mobiles  et  d’estimer  leurs  parametres  de  vitesse.  Le  systeme  radar  a  deux 
ouvertures  presente  toutefois  une  grave  limitation  en  matiere  de  GMTI,  en  raison  de  la 
contamination  des  signaux  de  cibles  mobiles  par  le  clutter.  Cette  limitation  a  stimule 
la  recherche  sur  les  systemes  multicanaux  utilisant  des  techniques  de  basculement  ou 
de  commutation  d’ouvertures,  grace  auxquelles  le  radar  a  deux  ouvertures  materielles 
se  comporte  comme  un  systeme  a  ouvertures  multiples  en  creant  des  centres  de  phase 
virtuels.  Le  present  document  technique  propose,  examine  et  evalue  trois  conceptions 
de  modes  multicanaux  pour  le  radar  d’imagerie  de  type  RADARSAT-2. 

Resultats: 

Techniques  multicanaux 

Trois  modes  a  centres  de  phase  multiples  sont  proposes  pour  MODEX  de  RADARS AT- 
2.  Le  premier  est  un  mode  d’emission  par  basculement  a  trois  canaux,  dans  lequel  le 
radar  emet  en  alternance  dans  les  sous-ouvertures  avant  et  arriere  et  reoit  les  echos  en 
utilisant  simultanement  les  deux  sous-ouvertures.  Les  deuxieme  et  troisieme  modes 
font  respectivement  appcl  a  des  techniques  de  reception  par  basculement  a  trois  ca¬ 
naux  et  a  quatre  canaux.  Dans  ces  modes,  le  radar  emet  dans  l’ouverture  complete, 
rnais  reoit  les  echos  dans  des  sous-sections  differentes  de  l’ouverture,  au  moyen  d’im- 
pulsions  en  alternance.  Le  mode  d’emission  par  basculement  a  trois  canaux  offre  la 
merne  ligne  de  base  longitudinale  quc  le  systeme  original  a  deux  ouvertures.  Toutefois, 
pour  les  modes  de  reception  par  basculement  a  trois  canaux  et  a  quatre  canaux,  la 
ligne  de  base  correspond  a  la  rnoitie  de  cclle  du  systeme  a  deux  ouvertures.  Line  ligne 
de  base  longitudinale  plus  courte  donne  habituellement  une  sensibilite  plus  basse  aux 
cibles  mobiles  lcntes.  Par  ailleurs,  un  degre  de  liberte  superieur  (nombre  de  canaux 
disponibles)  peut  ameliorer  le  rendement  de  detection  et  compenser  ainsi  la  reduction 
de  la  ligne  de  base. 
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Compromis  entre  la  largeur  du  faisceau  et  le  niveau  des  lobes 
lateraux 


Differentes  combinaisons  de  taille  des  ouvertures  d’emission  et  de  reception  produisent 
differents  diagrammes  d’antenne  bidirectionnels  et  des  variations  correspondantes 
touchant  le  niveau  des  lobes  lateraux  et  la  largeur  du  faisceau.  La  plus  grande  largeur 
de  faisceau  (mais  pas  le  niveau  le  plus  eleve  des  lobes  lateraux)  s’obtient  pour  une 
emission  sur  demi-antenne  et  une  reception  sur  demi-antenne  (c.-a-d.  dans  le  mode 
d’emission  par  basculement).  Le  diagramme  d’antenne  produit  pour  une  emission 
sur  antenne  complete  et  une  reception  sur  quart  d’antenne  (c.-a-d.  dans  lc  mode  de 
reception  par  basculement  a  quatre  canaux)  ne  differe  pas  beaucoup  de  remission  sur 
antenne  complete  et  de  la  reception  sur  demi-antenne  (c.-a-d.  le  mode  de  reception  a 
deux  ouvertures). 

Amelioration  de  I’estimation  de  la  phase  interferometrique 
longitudinale  et  de  la  detection  des  cibles  mobiles 

II  est  demontre  que  le  rendement  de  GMTI  multicanaux  fonde  sur  de  simples  mesures 
de  detection  avec  antennes  a  centre  de  phase  deplane  (DPCA)  s’ameliore  en  fonction 
du  nombre  de  canaux  disponibles  pour  un  espacement  constant  des  centres  de  phase. 
De  meme,  il  s’avere  que  l’estimation  de  la  phase  interferometrique  longitudinale,  qui 
est  proportionnelle  a  la  vitesse  radiate  d’une  cible  mobile,  augmente  de  precision 
lorsque  le  nombre  de  canaux  s’accroit.  Ces  resultats  s’appliquent  seulement  au  cas 
de  la  ligne  de  base  constante.  L ’augmentation  du  nombre  de  canaux,  au  detriment 
d’une  reduction  de  la  ligne  de  base,  peut  ne  pas  entrainer  d’augmentation  globale  du 
rendement  de  GMTI. 

Annulation  de  clutter  multicanaux 

La  production  de  trois  ou  quatre  centres  de  phase  dans  un  systeme  ROS  a  deux 
canaux  souleve  une  question  essentielle  :  ces  centres  de  phase  virtuels  produits  par 
basculement  d’ouverture  peuvent-ils  s’utiliser  efhcacement  pour  supprimer  le  clutter 
de  fond  fixe  ?  Des  experiences  de  simulation  imitant  les  trois  modes  multicanaux 
proposes  demontrent  clairement  la  possibilite  d’une  annulation  efhcace  du  clutter 
fixe  au  moyen  des  centres  de  phase  produits  par  le  basculement  de  sous-ouverture,  de 
faon  que  la  phase  interferometrique  d’une  cible  mobile  puisse  s’estimer  avec  plus  de 
precision. 

Portee:  La  possibilite  d’utiliser  la  technique  de  basculement  ou  de  commutation 
d’ouverture  pour  l’application  MODEX  multicanaux  de  RADARSAT-2  est  demontree 
ici  pour  la  premiere  fois  au  moyen  d’une  simulation.  Bien  que  l’application  de  cette 


VI 
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technique  a  la  GMTI  multicanaux  ne  soit  pas  nouvelle,  la  ”  preuve  ”  de  conception  en 
est  presentee  ici.  Ce  resultat  concluant  souleve  la  reelle  possibility  de  modes  MODEX 
multicanaux  pour  RADARSAT-2,  ce  qui  n’etait  pas  considere  comme  possible  durant 
les  etapes  initiales  de  conception  de  MODEX. 

Recherches  futures:  Les  recherches  futures  devraient  comprendre  des  experiences 
de  simulation  complete  au  moyen  de  SBRMTISIM  (simulateur  d’indication  de  cibles 
mobiles  de  radar  spatial)  en  vue  d’examiner  les  avantages  et  les  inconvenients  de 
chacun  des  modes  proposes. 
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Figure  1:  Concepts  of  increasing  the  spatial  diversity  for  RADARSAT-2, 
allowing  the  two-channel  SAR  system  to  effectively  operate  as  a 
three  or  four  channel  radar.  This  can  either  be  accomplished  by 
transmitter  toggling  between  subsequent  pulses  or  via  receiver 
excitation  and  toggling  schemes .  2 

Figure  2:  The  antenna  patterns  of  different  transmit  and  receive 

sub-aperture  combinations  showing  the  different  resulting 
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sub-apertures  (eight  sub-panels  each)  and  returns  are  received 
using  16  sub-panels.  The  return  signals  received  at  the  16 
sub-panels  are  added  in  such  a  way  as  to  form  three  sub-apertures 
with  their  respective  phase  centres .  7 

Figure  5:  Schematic  showing  the  simulation  of  a  multi-aperture  antenna: 
radar  signal  is  transmitted  at  full  aperture  and  returns  are 
received  using  16  sub-panels.  The  return  signals  received  at  the 
16  sub-panels  are  added  in  such  a  way  as  to  form  three 
sub-apertures  or  phase  centres .  8 

Figure  6:  Schematic  showing  the  simulation  of  a  multi-aperture  antenna: 

the  radar  signal  is  transmitted  using  the  full  aperture  and  returns 
are  received  using  16  sub-panels.  The  return  signals  received  at 
the  16  sub-panels  are  added  in  such  a  way  as  to  form  four  phase 
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1  Introduction 


Canada’s  RADARSAT-2  commercial  spaceborne  SAR,  to  be  launched  in  2007,  can 
be  operated  in  an  experimental  mode  (called  MODEX  for  Moving  Object  Detection 
Experiment)  which  will  allow  the  full  antenna  to  be  broken  into  two  sub-apertures 
with  two  parallel  receivers  to  define  two  independent  data  channels  [2],  These  two 
sub-apertures  in  the  flight  direction  can  be  used  to  record  two  echoes  (the  dual-receive 
mode),  one  from  each  wing  for  every  pulse  transmitted  using  the  full  antenna.  The 
two  apertures  make  it  possible  to  detect  targets  with  non-zero  across-track  velocity  by 
providing  two  identical  views  of  the  observed  scene  obtained  at  slightly  different  times. 
Recently,  several  new  detection  techniques  for  two-channel  SAR  ground  moving  target 
indication  (GMTI),  including  their  statistical  performance,  have  been  proposed,  e.g. 
[3].  It  is  well  known  that  quasi-optimum  detection  and  parameter  estimation  can  be 
achieved  if  at  least  three  physical  phase  centres  are  available  [4]  [5]  [6]:  Two  degrees- 
of-freedom  are  not  enough  to  simultaneously  suppress  the  clutter  and  to  accurately 
estimate  the  target’s  properties,  such  as  velocity  and  location. 

Although  the  design  of  RADARSAT-2  only  permits  the  simultaneous  use  of  two 
receivers,  an  alternating  transmit  mode  (or  the  toggle-transmit  mode  [10])  is  antici¬ 
pated  to  provide  a  third  independent  phase  centre  [11].  In  this  toggle-transmit  mode, 
pulses  are  alternately  transmitted  from  the  fore  and  aft  apertures  and  are  simulta¬ 
neously  received  using  both  apertures.  This  mode  not  only  allows  greater  separation 
of  the  two-way  phase  centres  in  the  along-track  direction  but  also  the  possibility  of 
generating  a  third  phase  centre  for  three-aperture  GMTI.  In  addition  to  this  mode  of 
operation,  the  very  flexible  programming  capabilities  of  the  RADARSAT-2  antenna 
allow  the  two-channel  SAR  system  to  operate  practically  like  a  three  or  four  channel 
radar  using  schemes  of  well-designed  transmitter /receiver  excitation. 

This  Technical  Memorandum  begins  by  describing  concepts  of  RADARSAT-2  multi¬ 
channel  modes,  followed  by  analysis  and  demonstration  of  the  performance  gains  for 
SAR-GMTI  by  the  exploitation  of  additional  channels.  The  results  of  simulation 
experiments  are  presented. 

2  Multi-Channel  Concepts 


Several  multi-aperture  concepts  have  been  explored.  These  are  illustrated  in  Fig.  1. 
By  doubling  the  pulse  repetition  frequency  (PRF)  of  the  radar,  it  is  possible  to 
generate  third  and  fourth  phase  centres  using  the  three  basic  schemes  as  illustrated. 
The  Erst  one,  Fig.  1(a),  uses  the  so-called  toggle-transmit  mode,  in  which  the  radar 
alternately  transmits  using  its  fore  and  aft  sub-apertures  and  receives  the  echoes  with 
both  sub-apertures  simultaneously.  This  scheme  allows  the  generation  of  three  phase 
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Figure  1 :  Concepts  of  increasing  the  spatial  diversity  for  RADARSAT-2,  allowing  the 
two-channel  SAR  system  to  effectively  operate  as  a  three  or  four  channel  radar.  This 
can  either  be  accomplished  by  transmitter  toggling  between  subsequent  pulses  or  via 
receiver  excitation  and  toggling  schemes. 


2 


DRDC  Ottawa  TM  2006-209 


centres  as  indicated  by  the  black  triangles  in  Fig.  1.  The  two  “side”  phase  centres  are 
used  to  sample  radar  echoes  at  half  the  rate  of  the  middle  phase  centre.  Therefore, 
in  order  to  avoid  the  aliasing  of  the  “clutter  band,”  it  is  necessary  to  double  the  pulse 
repetition  frequency  (PRF)  of  the  radar.  However,  the  PRF  is  not  the  only  parameter 
of  concern.  Since  the  radar  is  now  transmitting  at  half  the  antenna  width,  there  is  a 
broadening  of  the  physical  beamwidth,  which  increases  the  clutter  bandwidth  and,  in 
turn,  the  aliasing  of  the  clutter  band.  This  beam  broadening  will  be  the  subject  of  the 
next  section.  For  now,  it  is  noted  that  there  may  be  a  small  aliasing  with  this  mode  of 
operation  given  that  the  Doppler  spectrum  as  given  by  the  3-dB  beamwidth  is  about 
3752  Hz  whereas  the  maximum  RADARSAT-2  PRF  is  only  3800  Hz.  The  along- 
track  baseline  of  this  approach  is  twice  as  large  as  that  for  the  other  two  schemes. 
The  longer  baseline  will  permit  the  detection  of  slower  moving  targets,  because  the 
along-track  interferometric  phase  of  a  moving  target  is  proportional  to  the  baseline. 

The  second  and  the  third  schemes,  instead  of  using  different  sub-apertures  on  trans¬ 
mit,  use  different  sub-apertures  on  receive.  This  is  accomplished  by  “turning  off’ 
different  regions  of  the  antenna  on  receive,  so  that  echoes  are  sampled  by  different  sub¬ 
panels  of  the  antenna  in  an  alternating  manner,  as  illustrated  in  Figs.  l(b)&(c).  The 
second  scheme  generates  three  phase  centres:  The  middle  phase  centre  is  generated 
by  summing  the  two  half  portions  of  the  fore  and  aft  sub-apertures  of  RADARSAT-2. 
The  along-track  baseline  (the  distance  between  phase  centres)  of  this  approach  is  half 
of  that  in  the  two-aperture  along-track  interferometry  (ATI)  but  offers  greater  spatial 
diversity  than  the  two-aperture  system.  It  provides  the  advantage  of  an  extra  degree 
of  freedom  for  clutter  cancellation  resulting  in  a  potential  improvement  in  target  pa¬ 
rameter  estimation.  Again,  the  PRF  must  be  doubled  to  avoid  aliasing.  The  third 
approach  to  multi-channel  GMTI,  Fig.  1(c),  is  very  similar  to  the  second  scheme  but 
produces  four  phase  centres  with  the  same  along-track  baseline  as  the  second  scheme. 
Different  regions  of  sub-apertures  are  excited  alternately  at  receive  such  that  a  pair 
of  phase  centres  are  generated  for  each  alternating  pulse.  The  PRF  is  doubled  to 
minimize  the  aliasing  of  the  clutter  band.  It  should  also  be  noted  that  there  is  a 
price  to  pay  for  increasing  the  PRF,  since  this  reduces  the  width  of  the  swath  in  slant 
range.  This  could  impact  the  utility  of  this  scheme  for  different  applications. 


3  Beamwidth  and  Sidelobe  Level  Trade-Offs 


In  this  section  the  trade-offs  in  beamwidth  and  sidelobe  level  for  the  three  proposed 
multi-aperture  modes  are  examined.  When  partitioning  a  single  aperture  into  smaller 
parts,  the  two-way  beamwidth  (or  antenna  pattern)  changes.  This  can  be  easily 
seen  in  Fig.  2.  The  beam  patterns  are  computed  using  the  MATLAB  code  ”beam- 
Forming2.m”  (given  in  Appendix  C).  The  different  combinations  of  transmit /receive 
aperture  sizes  yield  different  two-way  antenna  patterns,  including  increased  sidelobe 
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Figure  2:  The  antenna  patterns  of  different  transmit  and  receive  sub-aperture  com¬ 
binations  showing  the  different  resulting  beamwidths  and  sidelobe  levels.  “F-F”:  full 
aperture  transmit  and  full  aperture  receive,  “F-H”:  full  aperture  transmit  and  half 
aperture  receive,  “F-Q”:  full  aperture  transmit  and  quarter  aperture  receive,  and 
“ H-H” :  half  aperture  transmit  and  half  aperture  receive. 
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levels  and/or  increased  beamwidths.  The  narrowest  beamwidth  and  lowest  sidelobe 
levels  are  obtained,  as  expected  ,  using  the  full  aperture  on  transmit  and  on  receive. 
The  widest  beamwidth  (but  not  the  highest  sidelobe  levels)  is  generated  using  half 
of  the  antenna  on  transmit  and  half  of  the  antenna  on  receive.  It  is  interesting  to 
note  that  the  antenna  patterns  generated  from  the  ”F-H”  (full  aperture  transmit  and 
half  aperture  receive)  and  the  ”F-Q”  (full  aperture  transmit  and  quarter  aperture 
receive)  do  not  differ  significantly  in  the  main  beam;  only  their  sidelobe  levels  are 
different  (about  3  dB  difference).  This  is  good  news  since  there  is  not  a  significant 
loss  in  buying  these  extra  degrees  of  freedom. 

The  3-dB  main-beam  induced  clutter  Doppler  spread  can  be  computed  from  the  radar 
platform  velocity  va  projected  onto  the  radar-to-clutter  radial  direction.  Assuming 
zero-angle  beam  squint,  this  is  given  by  [1] 

A/d  =  2/d  =  2(2|)=4^,  (1) 

where  vs  and  vr  are  the  spacecraft  velocity  and  the  radial  component  of  this  velocity, 
respectively,  2 ip  is  the  3-dB  beamwidth  angle,  and  A  is  the  wavelength. 

4  Moving  Target  Signal  Model 


Simulation  experiments  have  been  carried  out  to  validate  the  multi-aperture  concept 
for  the  RADARSAT-2  MODEX.  For  simplicity,  the  flat-earth  or  airborne  geometry 
is  chosen  for  the  simulation,  since  the  “proof’  of  the  concept  is  independent  of  the 
imaging  geometry.  The  phase  history  of  a  moving  target  in  a  flat-earth  geometry,  as 
shown  in  Fig.  3,  can  be  modelled  as 

<p(t)  =  —2kR{t) 

I -  (2) 

=  —2k\J[xo  +  (vx0  -  va)t  +  ax0t2/ 2]2  +  [yQ  +  vy0t  +  ay0t2/2}2  +  h 2, 

where  R(t)  is  the  range  history  of  the  moving  target,  k  =  2n/\  (A  is  wavelength),  and 
v  and  a  are  target  velocity  and  acceleration,  respectively.  Subscript  “0”  denotes  values 
at  t  =  0,  and  x  and  y  are  the  along-track  and  across-track  components,  respectively. 
The  target  velocity  and  position  at  broadside  time  tb  are  therefore 

Vxb  ArO  T  Ojxq  tb, 

Vyb  VyO  T  ttyotbi  ^2^ 

Xb  Vatb  2^0  T  ^xO tb  T 
Ub  —  y  o  +  Vyo  tb  +  ayotl/2. 

The  second  last  equation  of  (3)  is  a  result  of  the  definition  of  “broadside,”  where  the 
azimuth  position  of  the  radar  coincides  with  that  of  the  moving  target. 
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Figure  3:  A  target  moving  with  constant  acceleration  in  a  Hat-earth  geometry  with 
components  of  motion  in  the  x  and  y  directions.  The  radar  moves  in  the  positive  x 
direction. 


Assuming  a  constant-acceleration  target  (i.e.  constant  ax o  and  ay o),  the  Taylor  series 
expansion  of  (2)  about  tb  can  be  written  as 


<p(t)  ft  -2 k{Rb  +  Vyb'fit  -  tb) 


+  2R^ei  +  ybayo)(t  —  tb)2 

T  _  \(Xxb  Va)cLx 0  T  Xyb(lyQ ( 1  y  )]  (t  tfe)  +...}, 

ZKb 

where  is  the  slant  range  to  the  target  at  broadside  and 


(4) 


7  =  Vb/Rbr 

2  /  \2  ,  2  /-i  2\ 

Vrel  —  \vxb  Va)  T  VybA  7  )• 

The  detailed  derivation  of  (4)  is  given  in  Appendix  B.  The  “slow-time”  moving  target 
signal  can,  therefore,  be  modelled  as 


s  =  rect 


t-tb 

T 


(6) 


which  is  a  finite-time  linear  chirp  with  duration  T  centred  at  time  t  =  tb-  For  actual 
signals,  the  “rect”  function  would  be  modulated  by  the  antenna  pattern. 
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Figure  4:  Schematic  showing  the  simulation  of  a  multi-aperture  antenna:  the  radar 
pulse  is  transmitted  at  alternating  fore  and  aft  sub-apertures  (eight  sub-panels  each) 
and  returns  are  received  using  16  sub-panels.  The  return  signals  received  at  the 
16  sub-panels  are  added  in  such  a  way  as  to  form  three  sub-apertures  with  their 
respective  phase  centres. 

5  Experiments 

5.1  Simulation  Design 

In  order  to  test  and  validate  the  multi-channel  concepts  obtained  by  sub-aperture 
toggling  (or  switching)  for  a  physical  two-channel  radar,  simulation  experiments  can 
be  designed  in  MATLAB.  These  mimic  what  actually  takes  place  in  the  sampling  of 
radar  echoes  using  the  temporal  diversity  strategies  described  in  the  previous  section 
and  the  subsequent  GMTI  signal  processing. 

The  first  simulation  considers  the  toggle-transmit  mode  of  the  RADARSAT-2  MODEX. 
The  simulator  architecture  is  illustrated  in  Fig.  4.  The  three  phase  centres  are  gen¬ 
erated  by  alternately  transmitting  radar  pulses  using  the  fore  and  aft  apertures  (as 
shown  in  solid  and  dash-lined  arrows)  and  sampling  the  echoes  using  sixteen  sub¬ 
panels  (shown  as  green  rectangles),  as  for  the  sixteen  sub-panels  on  RADARSAT-2. 
The  signals  received  by  the  sixteen  sub-panels  are  then  coherently  added  to  form 
channels  1  and  2  for  the  odd  pulses  and  channels  2  and  3  for  the  even  pulses.  The 
positions  of  the  actual  phase  centres  on  the  antenna  are  shown  as  red  dots  in  Fig. 
4.  The  MATLAB  code  simulator  ”airMultiChan_toggle.m”  is  used  to  validate  the 
concept  and  is  given  in  Appendix  E. 
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Figure  5:  Schematic  showing  the  simulation  of  a  multi-aperture  antenna:  radar  signal 
is  transmitted  at  full  aperture  and  returns  are  received  using  16  sub-panels.  The 
return  signals  received  at  the  16  sub-panels  are  added  in  such  a  way  as  to  form  three 
sub-apertures  or  phase  centres. 

The  second  multi-channel  simulation  architecture  is  illustrated  in  Fig.  5.  The  radar 
signal  is  transmitted  using  the  full  aperture  and  the  echoes  are  sampled  using  16 
sub-panels.  The  signals  from  the  16  sub-panels  are  added  in  such  a  way  so  as  to 
form  three  phase  centres  (shown  as  red  dots  in  the  figure).  This  architecture  allows 
the  simulation  of  a  true  multi-channel  radar  as  well  as  the  emulation  of  a  multi¬ 
channel  radar  via  temporal  diversity  (or  switching).  The  three  phase  centres  can  be 
generated  by  sampling  the  echoes  using  three  different  sub-sections  of  the  antenna  as 
shown  in  the  figure.  Since  RADARSAT-2  only  has  two  channels  for  the  sampling  of 
return  echoes,  the  apertures  for  phase  centres  1  and  3  must  be  sampled  for  the  odd 
numbered  pulses  and  those  for  phase  centre  2  for  the  even  numbered  pulses.  The 
PRF  must  be  doubled  in  order  to  maintain  the  original  sampling  rate  (for  each  phase 
centre)  to  avoid  aliasing. 

The  third  multi-channel  simulation  architecture  (illustrated  in  Fig.  6)  is  very  similar 
to  the  previous  one.  In  the  same  way,  the  radar  pulse  is  transmitted  at  full  aperture 
and  return  echoes  are  sampled  using  16  sub-panels.  Then  the  signals  from  the  16 
sub-panels  are  added  coherently  to  form  four  distinct  phase  centres  as  shown.  The 
simulator  permits  the  simulation  of  a  true  four-channel  radar  as  well  as  the  emulation 
of  a  four-channel  radar  via  sub-aperture  switching.  With  only  two  physical  receiver 
channels  of  RADARSAT-2,  it  is  necessary  to  sample  phase  centres  1  and  3  on  the 
odd  numbered  pulses  and  phase  centres  2  and  4  on  the  even  numbered  pulses.  Again, 
the  PRF  must  be  at  least  doubled  in  order  to  keep  the  original  sampling  rate  at  each 
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Figure  6:  Schematic  showing  the  simulation  of  a  multi-aperture  antenna:  the  radar 
signal  is  transmitted  using  the  full  aperture  and  returns  are  received  using  16  sub- 
panels.  The  return  signals  received  at  the  16  sub-panels  are  added  in  such  a  way  as 
to  form  four  phase  centres. 

phase  centre.  The  doubling  of  the  PRF  may  not  be  sufficient  due  to  the  changes  in 
the  antenna  patterns  in  the  toggle  modes. 


5.2  Toggle-Transmit  Mode 


RADARSAT-2  MODEX,  as  mentioned  earlier,  is  capable  of  emulating  a  third  phase 
centre  by  operating  in  the  togglc-transmit  mode,  thereby  allowing  the  three-channel 
processing  of  ground  moving  targets.  With  alternating  transmit  and  simultaneous 
receive  using  fore  and  aft  apertures,  the  effective  two-way  3  dB  beamwidth  is  about 
0.4°  in  azimuth  resulting  in  a  clutter  bandwidth  of  1876  kHz.  As  the  toggle-transmit 
mode  requires  a  PRF  that  is  at  least  twice  the  clutter  bandwidth,  the  maximum  PRF 
for  RADARSAT-2  of  3.8  kHz  is  only  marginally  above  the  Nyquist  frequency.  With 
the  third  phase  centre,  the  interfering  clutter  can  be  cancelled  before  forming  the 
interferogram  by  first  subtracting  (via  DPCA)  the  signals  in  channel  2  from  those 
in  channel  1  and  those  in  channel  3  from  those  in  channel  2,  and  then  forming  the 
interferogram  from  the  two  resulting  clutter-cancelled  “channels.”  For  channels  1  and 
2,  the  phase  of  the  DPCA  output  signal  can  be  shown  to  be 


Vh2 


2n(Ri  +  i?2 )  .  7r 

A  2’ 


(7) 


where  R±  and  i?2  are  the  distances  from  the  phase  centres  1  and  2  to  the  target  of 
interest,  respectively.  Similarly  for  channels  2  and  3,  1^23  =  27r(i?2  +  Rs)/^  ±  7t/2. 
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The  resulting  interferometric  phase  is  thus 


2  7r(i?i  —  R3) 

A 


47r(i?i  -  R2 ) 

A 


47rtyr 

A 


(8) 


The  phase  of  the  three-channel  ATI,  in  the  toggle-transmit  mode,  is  exactly  the 
same  as  that  for  the  two-channel  ATI,  except  the  clutter  has  now  been  cancelled  via 
the  DPCA  subtraction.  All  the  parameter  estimation  algorithms  previously  derived 
for  the  two-channel  ATI  can  now  be  applied  in  exactly  the  same  manner  without 
any  modification.  A  simulation  for  a  moving  target  with  vx  =  —26.63  m/s  and 
vr  =  —10.79  m/s  in  a  background  of  Gaussian  clutter  (and  using  a  matched  filter 
perfectly  matched  to  the  target  motion)  shows  that  the  target’s  ATI  phase  deviates 
strongly  from  the  ideal  (uncorrupted)  phase  value  of  70.5°  when  the  two-phase-centre 
ATI  approach  is  used  (see  the  blue  curves  in  Fig.  7).  This  is  believed  to  be  due  to 
corruption  of  the  target  signal  by  the  overlapping  stationary  clutter.  The  ATI  phase 
bias  is  a  function  of  the  signal-to-clutter  ratio  (SCR).  With  the  three-phase-centre 
ATI,  on  the  other  hand,  the  clutter-cancelled  target  signals  yield  an  unbiased  ATI 
phase  (as  shown  in  the  red  curves  of  Fig.  7). 


5.2.1  Can  the  Toggle-Transmit  Mode  Really  Produce  Three  Phase 
Centres? 

Although  generating  three  or  four  phase  centres  from  a  two-channel  SAR  system  may 
seem  at  first  magical,  German  researchers  have  been  using  this  aperture  switching 
concept  to  increase  spatial  diversity  of  a  two-channel  radar  system  for  many  years 
[7]  [8]  [9].  The  key  question  that  needs  to  be  answered  is,  “Does  the  toggle-transmit 
mode,  derived  from  the  two  physical  receiver  channels  of  RADARSAT-2,  really  pro¬ 
vide  an  independent  third  phase  centre  useable  for  clutter  suppression?” 

The  simulator  “airMultiChan_toggle.m”  is  used  to  generate  the  moving  target  signal 
[yx  =  3  and  vy  =  —4)  embedded  in  a  competing  clutter  signal.  It  then  samples  radar 
echoes  using  the  strategies  (to  be  used  in  the  actual  RADARSAT-2  signal  processor), 
which  are  described  earlier  and  generates  three  signal  channels  sampled  at  twice  the 
original  PRF.  The  clutter  signals  are  cancelled  using  the  procedure  described  above, 
followed  by  the  formation  of  the  interferogram  from  the  “clutter-free”  signals.  The 
results  of  the  simulator  output  are  shown  in  Figs.  8(a)  to  9(e). 

Fig.  8(a)  depicts  the  antenna  “pattern”  obtained  by  summing  the  signals  received 
by  eight  sub-panels  of  the  sixteen-panel  antenna.  Fig.  8(b)  shows  the  compressed 
moving  target  signal  embedded  in  a  clutter  background.  The  signal  compression 
is  carried  out  using  the  fractional  Fourier  transform  [18].  Figs.  8(c)  through  9(a) 
are  time-frequency  representations  (via  the  short-time  Fourier  transform)  of  the  fore 
aperture  signals,  the  fore  channel  minus  the  centre  channel  1,  the  centre  channel 
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Figure  7:  A  simulated  deterministic  target  moving  with  vx  =  —26.63  m/s  and  vr  = 
—  10.79  m/s  in  a  background  Gaussian  clutter.  The  angle  of  incidence  is  50°,  the 
matched  filter  is  matched  to  the  target  motion,  and  the  ideal  ATI  phase  is  70.5°.  The 
two-phase-centre  and  three-phase-centre  results  are  denoted  by  blue  and  red  lines, 
respectively.  The  SCR  is  calculated  before  azimuth  compression;  both  signal  and 
matched-filter  lengths  are  equal  to  the  4  dB  beam  width;  scene  coherence  p  =  0.95. 
(a)  SCR  =  -24.4  dB;  (b)  SCR  =  -21.6  dB;  (c)  SCR  =  -17  dB;  (d)  SCR  =  -6.0 
dB. 
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1  minus  the  aft  channel,  the  centre  channel  2  minus  the  aft  channel,  and  the  fore 
channel  minus  the  aft  channel.  As  can  be  seen,  the  simulation  results  clearly  show 
that  the  generation  of  the  third  independent  phase  centre  is  indeed  feasible  using 
the  togglc-transmit  scheme.  Figs.  8(c)-8(f)  clearly  demonstrate  that  the  fore  and  aft 
channel  clutter  signals  can  be  suppressed  using  the  third  (virtual)  phase  centre  (or 
channel  2).  The  interferograms  of  the  fore  and  aft  channels  are  initially  contaminated 
by  the  interfering  clutter  as  seen  in  Figs.  9(b)-9(c).  After  the  clutter  cancellation,  the 
interferogram  is  well  defined  and  its  phase  agrees  with  the  true  phase  of  the  moving 
target. 


5.3  Toggle-Receive  Modes 

Instead  of  changing  the  antenna  pattern  on  transmission  (i.e.  via  excitation  of  differ¬ 
ent  parts  of  the  antenna),  the  alternative  is  to  vary  the  regions  of  antenna  excitation 
at  reception.  Radar  pulses  are  transmitted  using  the  full  antenna  aperture,  and  re¬ 
turn  echoes  are  received  with  sub-sections  of  the  antenna  in  order  to  provide  greater 
spatial  diversity  for  more  efficient  clutter  cancellation. 

5.3.1  Improvement  in  ATI  Phase  Estimation  and  Detection 

The  key  question  to  be  answered  is  how  these  different  modes  of  antenna  excitation 
(with  an  increased  number  of  available  independent  channels)  affect  or  improve  the 
measurement  of  a  moving  target’s  interferogram?  In  the  following,  simulations  of  the 
moving  target’s  interferometric  phase  (contaminated  by  background  clutter)  for  two, 
three,  and  four  channel  SAR  systems  are  performed  and  ATI  phase  probability  density 
functions  (PDFs)  as  a  function  of  the  number  of  channels  (showing  improvement  on 
accuracy  of  phase  measurements)  are  shown  in  Fig.  10.  As  before,  the  ATI  phase  of 
the  stationary  clutter  is  centred  around  the  zero  value  and  the  phase  PDFs  of  a  target 
moving  at  vx  =  3  m/s  and  vy  =  4  m/s  are  also  shown.  The  conventional  two-channel 
system,  computed  using  only  channels  1  and  2  or  channels  2  and  3,  yields  biased  ATI 
phase  measurements  (labelled  as  ‘12  &  23’  on  the  figure).  The  ATI  phase  measured 
from  the  three-channel  SAR,  computed  by  forming  an  interferogram  between  the 
DPCA  of  channels  1  and  2  and  the  DPCA  of  channels  2  and  3  or  between  the  DPCA 
of  channels  2  and  3  and  the  DPCA  of  channels  3  and  4  (labelled  as  ‘123  &  234’), 
shows  no  phase  biases.  However,  the  measured  phase  variance  can  be  quite  large 
depending  on  the  degree  of  clutter  cancellation  or  the  clutter  coherence  between  the 
channels.  When  all  four  channels  are  used,  the  ATI  phase  measurements  obtained 
by  summing  the  interferogram  measurements  from  channels  1  to  3  and  channels  2  to 
4  give  an  unbiased  ATI  phase  with  a  noticeably  smaller  variance  (labelled  as  ‘ATI 
sum’)).  Also  shown  are  the  ATI  phase  measurements  computed  using  the  following 
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Figure  8:  Three-channel  GMTI  using  toggle-transmit  mode:  alternating  transmit  at 
fore  and  aft  half  sub-apertures  and  receive  simultaneously  with  two  half  sub-apertures. 
The  fore  and  aft  channels  are  sampled  at  half  the  rate  of  the  centre  channel.  To  avoid 
aliasing,  the  PRF  is  doubled,  (a)  Normalized  antenna  patterns,  (b)  fore  channel 
signals  (target  plus  clutter)  after  fractional  Fourier  transform  compression,  time- 
frequency  plots  of  (c)  fore  aperture  signals,  (d)  fore  channel  minus  centre  channel  1, 
(e)  centre  channel  1  minus  aft  channel,  and  (f)  centre  channel  2  minus  aft  channel. 
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Measured  ATI  Phase  1  &  2  =  -88.2099 


Measured  ATI  Phase  2  &  3  =  -100.4629 
90  80000 


Measured  ATI  Phase  dpcal  &  dpca2  (center  1)  =  -85.6752 
90  150000 


(c) 


(d) 


Measured  ATI  Phase  dpcal  &  dpca2  (center  2)  =  -85.4528 
90  150000 


Theoretical  ATI  Phase  =  -85.3573 

(e) 


Figure  9:  The  same  mode  as  in  the  previous  figure,  (a)  Time- frequency  plots  of 
fore  channel  minus  aft  channel  and  polar  plots  of  interferogram  between  (b)  fore 
and  centre  1  channels,  (c)  centre  1  and  aft  channels,  (d)  DPCA  of  fore  and  centre  1 
channels  and  DPCA  of  centre  1  and  aft  channels,  and  (e)  DPCA  of  fore  and  centre  2 
channels  and  DPCA  of  centre  2  and  aft  channels. 
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Figure  10:  ATI  phase  PDFs  for  two,  three,  and  four  channel  systems. 


steps: 


DPCA123  =  DPCA12  -  DPCA23, 

DPCA234  =  DPCA23  -  DPCA3A,  (9) 

ATI  =  ( DPCA123){DPCA23A )*, 

where  DPCA12,  DPCA23 ,  and  DPCA3A  are  output  signals  obtained  by  subtracting 
signals  between  channels  1  and  2,  channels  2  and  3,  and  channels  3  and  4,  respectively. 
As  can  been  in  the  figure  (labelled  as  ‘2  DPCAs’),  the  ATI  phase  computed  in  this 
manner  shows  a  slight  bias  in  the  phase  measurement  and  a  significantly  larger  vari¬ 
ance.  It  is  clear  that  this  is  not  the  way  to  compute  the  ATI  phase.  There  are  other 
ways  of  computing  the  ATI  phase  using  the  four  channels  available  to  us.  However, 
they  are  not  explored  further  in  this  Technical  Memorandum  (further  information  is 
included  in  “spacePDFgenerate_dpca.m”  in  Appendix  D). 

The  accuracy  of  ATI  phase  measurements  directly  impacts  estimation  of  the  target 
parameters.  But  how  does  increasing  the  number  of  channels  improve  the  probability 
of  detection?  This  question  can  be  examined  from  the  perspective  of  simple  DPCA 
clutter  cancellation.  In  the  following,  we  compare  the  signal-to-clutter  ratio  PDFs  for 
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DPCA  SCR  (dB) 


Figure  11:  DPCA  signal  magnitudes  for  two,  three,  and  four  channel  systems. 


two,  three  and  four  channel  SAR.  It  is  possible  to  cancel  the  clutter  via  DPCA  using 
various  combinations  of  DPCA  substruction  afforded  by  the  multiple  channels.  The 
first  is  the  conventional  two-channel  approach,  which  is  to  subtract  one  channel  from 
the  other  co-registered  channel  to  cancel  clutter.  The  second  uses  three  channels  to 
form  two  DPCAs:  one  from  channels  1  and  2  and  the  other  from  channel  2  and  3. 
Then  the  two  DPCA  outputs  are  coherently  added.  The  third  approach  uses  all  four 
channels  to  form  DPCAs  by  subtracting  channels  1  and  2,  channels  2  and  3,  and 
channels  3  and  4.  Then  all  three  DPCA  outputs  are  coherently  added.  Similarly, 
the  fourth  approach  uses  all  four  channels,  but  instead  of  forming  DPCA  using  the 
neighboring  channels,  the  DPCAs  are  computed  by  subtracting  channels  1  and  3  and 
channels  2  and  4.  These  are  then  coherently  added  together  to  form  a  new  DPCA 
output.  The  resulting  SCR  PDFs  for  the  above  four  approaches  are  shown  in  Fig.  11. 
As  expected,  the  SCR  increases  with  increasing  number  of  channels.  It  is  interesting 
to  observe  that  there  is  a  slight  variation  between  SCR  PDFs  for  the  third  and  fourth 
approaches.  It  is  not  clear  at  the  present  time  why  this  is  so. 
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5.3.2  Three-Channel  Toggle-Receive  Mode 

A  moving  target  embedded  in  background  stationary  clutter  is  simulated  for  the  three- 
channel  togglc-receive  mode,  as  illustrated  in  Fig.  5.  Shown  in  Figs.  12(b)  through 
13(a)  are  time-frequency  representations  of  signals  from  the  fore  channel,  the  fore 
channel  minus  the  centre  channel,  the  centre  channel  minus  the  aft  channel,  and 
the  fore  channel  minus  the  aft  channel.  Clearly,  clutter  cancellation  via  the  virtual 
channel  DPCA  subtraction  is  demonstrated  in  this  three-channel  togglc-receive  mode. 
Polar  plots  of  interferograms  between  the  fore  and  centre  channels,  the  centre  and 
aft  channels,  and  the  DPCA  of  the  fore  and  centre  channels  and  the  DPCA  of  the 
centre  and  aft  channels  are  also  shown  in  Figs.  13(b)  through  13(d).  Interferograms 
derived  from  channels  without  prior  clutter  cancellation  show  phase  contamination 
as  in  Figs.  13(b)  and  13(c).  The  virtual  centre  channel  is  again  demonstrated  to 
be  useful  in  clutter  suppression  as  seen  in  Fig.  13(d)  in  which  the  interferogram  is 
cleared  of  clutter  contaminants. 

5.3.3  Four-Channel  Toggle-Receive  Mode 

Simulation  results  for  the  four-channel  toggle-receive  mode  are  shown  in  Figs.  14(a) 
through  14(f).  They  are  similar  to  the  three-channel  toggle-transmit  and  togglc- 
receive  modes.  The  results  show  no  surprises.  Virtual  channels  generated  from  sub- 
aperture  switching  or  toggling  are  again  shown  to  be  effective  in  clutter  suppression  as 
seen  in  Figs  14(c)  through  14(f)  and  do  indeed  improve  the  interferogram  estimation 
accuracy  as  shown  in  Figs.  15(a)  through  15(e).  The  number  of  degrees  of  freedom  or 
spatial  diversity  required  to  carry  out  true  multi-channel  processing  of  moving  targets 
is  indeed  increased  via  the  proposed  sub-aperture  toggling  schemes. 


6  Conclusions 


This  Technical  Memorandum  proposes  and  describes  three  multi-channel  GMTI  modes 
for  the  RADARSAT-2  MODEX.  The  multi-channel  concepts  make  use  of  the  so-called 
aperture  switching  or  toggling  technique,  which  generates  independent  phase  centres 
and  increases  the  spatial  diversity  of  a  two-physical-channel  radar  to  become  a  three 
or  four  virtual-channel  radar.  Simulation  experiments  are  conducted  to  demonstrate 
the  feasibility  of  such  a  multi-channel  radar  system.  Results  clearly  show  that  virtual 
channels  can  be  exploited  to  improve  signal-to-clutter  ratio  (and,  therefore,  the  detec¬ 
tion  performance),  suppress  the  interfering  clutter,  and  enhance  the  target  parameter 
estimation  accuracy.  Other  considerations,  however,  need  to  be  taken  into  account 
when  applying  the  technique  such  as  a  reduced  range  swath,  beam  broadening,  and 
clutter  band  aliasing. 
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Figure  12:  Three-channel  GMTI  in  toggle-receive  mode  with  a  doubling  of  the  PRF. 
The  fore  aperture  consists  of  first  eight  sub-panels  of  the  antenna,  the  centre  aper¬ 
ture  consists  of  the  centre  eight  sub-panels,  and  the  aft  aperture  consists  of  the  last 
eight  sub-panels,  (a)  Normalized  antenna  patterns.  Time- frequency  plots  of  (b)  fore 
aperture  signal,  (c)  fore  channel  minus  centre  channel,  and  (d)  centre  channel  minus 
aft  channel. 
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Measured  ATI  Phase  1  &  2  =  -50.31 1 1 


Measured  ATI  Phase  2  &  3  =  -51.1556 


90  60000 


(c) 


90  60000 


(b) 


Measured  ATI  Phase  dpcal  &  dpca2  =  -42.9289 
90  50000 


(d) 


Figure  13:  The  same  mode  as  the  previous  figure,  (a)  Time- frequency  plot  of  fore 
channel  minus  aft  channel.  Polar  plots  of  the  interferogram  between  (b)  fore  and 
centre  channels,  (c)  centre  and  aft  channels,  and  (d)  DPCA  of  fore  and  centre  channels 
and  DPCA  of  centre  and  aft  channels. 
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(d) 
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(e) 


(f) 


Figure  14:  Four-channel  GMTI  in  toggle-receive  mode  with  a  doubling  of  the  PRF: 
transmit  at  full  aperture  and  receive  at  first  and  third  quarter  sub-apertures  on  odd 
pulses  and  receive  at  second  and  fourth  quarter  sub-apertures  on  even  pulses,  (a) 
Channel  1  signal  magnitude  (target  plus  clutter )  after  fractional  Fourier  transform 
compression,  time-frequency  plots  of  (b)  channel  1  signal,  (c)  channel  1  minus  channel 
2,  (d)  channel  3  minus  channel  4,  and  (f)  channel  1  minus  channel  3. 
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Measured  ATI  Phase  1  &  2  =  -42.0682 
90  20000 


Measured  ATI  Phase  2  &  3  =  -48.4552 
90  15000 


(a) 


(b) 


Measured  ATI  Phase  dpcal  &  dpca2  =  -42.8253 
90  10000 


(c) 


Measured  ATI  Phase  dpca2  &  dpca3  =  -42.6619 
90  15000 


(d) 


Measured  ATI  Phase  of  atil  +  ati2  =  -42.7433 
90  30000 


Theoretical  ATI  Phase  =  -42.6787 

(e) 


Figure  15:  The  same  mode  as  Fig.  15.  Polar  plots  of  the  phase  and  amplitude  of 
the  interferogram  between  (a)  channels  1  and  2,  (b)  channels  2  and  3,  (c)  DPCA  of 
channels  1  and  2  and  DPCA  of  channels  2  and  3,  (d)  DPCA  of  channels  2  and  3  and 
DPCA  of  channels  3  and  4,  and  (e)  by  adding  (c)  and  (d). 
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Annex  A:  Platform-Motion-Induced  Clutter  Doppler  Spread 

There  has  been  discussion  in  DRDC-Ottawa  Space-Based  Radar  Group  about  the  correct  way  to  compute  the  clutter 
Doppler  bandwidth  for  the  space-based  radar  geometry.  It  was  argued  that  a  radar’s  beam  sweeps  the  ground  at  slightly 
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Figure  A. 1 :  Spaceborne  geometry  for  clutter  bandwidth  calculation. 


r2  =  R2e  +  (h  +  Re)2  —  2Re(h  +  Re)  cos  9.  (A. 2) 

Taking  the  derivative  of  both  sides  of  (A. 2)  with  respect  to  time  t  yields 
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Annex  B:  Taylor  Series  Expansion  of  Radar-to-Target  Range 
History 


CO 

<D 

cd 

T— 1 

c/T 

’C 

CD 

_bb 

ed 

faO 

£ 

CD 

CO 

S-l 

fa 

a3 

X 

O 

.2 

a 

g 

£ 

£ 

i> 

P 

<P 

o 

CD 

hC 

0 

•a 

o' 

P 

.2 

CO 

kC 

II 

faO 

p 

Jh 

fao 

-+o 

P 

fa-1 

P 

CD 

c 

<P 

r2 

CD 

P 

"a; 

o 

CO 

CD 

P 

CD 

CD 

CD 

g 

hJO 

P 

'§ 

P 

fai 

> 

fan 

3 

f-H 

Oh 

Cl 

<P 

CD 

<M 

CO 

CD 

faO 

rT3 

QJ 

M 4 

<P 

fao 

<P 

+ 

o 

X 

.2 

CO 

X 

qp 

CD 

<N 

cf 

CD 

"Td 

P 

CD 

fan 

.2 

h£ 

fai 

CN 

-+o 

O 

CO 

CD  CO  . 
co  c  s 

-2  3  .2 

P  co 
co  ^  -r-j 

CD  P  § 

£  ®  a 
^  ^  -0 

£ 

-g  ;>  'o 
qj  +3  _2 
CL,  o  o 
So  CU  > 

OJ  (Oh 

"  S? 
a  ^  £? 
.2  «T  a 

43  +i  ^ 

Cd  HH  (V) 
S-H  QJ  p 

o  o  “ 
O  CO  . 


+  + 


e  e 
+  + 


b M  § 

CO  CD  CD 

p  rp  ^ 


S-S  I 

>>S  8 

‘o  ^  ^ 
X  co 
O  to 

CD  9  9 


w  .  'nU 

i?1  I 

P  P  P 

faO  CO 

CD  o  P 


co  .O  > 
co  eg  *C 

rt  *  | 

H  ®  !> 


„_,  bC  . 
rj  . 

3  o  2 


®  ^ 
CP  04  bC 

• — I  .+_*  H 

a  cu  o 

£  2  33 


DRDC  Ottawa  TM  2006-209 


27 


From  the  second  last  equation  of  (B.2),  —  x0  =  (vx0  -  va)tb  +  ax0tl/2.  Substituting  ‘—x0’  and  ‘yb  into  and  evaluating  (B.3) 
at  t  =  4  yields 

f>/i  \  VbiVyO  Q"yotb)  /-q  a  \ 

R(h)  = - 5 - =  ivyb,  (B.4) 
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where  v2rel  =  (vxb  -  vaf  +  v2b(l  -  72)  and  7  =  yb/Ri 
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Annex  C:  MATLAB  Code  ‘beamForming2.m’ 
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figure;  plot (theta*180/pi , 10*logl0 (abs (Et . *Er) )  -  EE_maxdB) 
xlabel (’ Angle  (degree)’);  ylabelC’ Amplitude  (dB) ’ ) 
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-32.3766;  ‘/calculated  from  orbital  paramters 
269.5196;  "/calculated  from  orbital  parameters 


delta_xO  =  -100;  °/0this  number  should  NOT  be  too  large  or  else  approximation  breaks  down 
tb  =  delta_xO/(Rs*ws-vx) ;  °/0this  is  a  negligible  approximation 
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-  (Rs*cos(phi)*(l  -  cos(ws*t))  +  2*d*sin(phi_s) *cos (ws*t) ) ) . ~2) ; 
s3u  =  b. *exp(-j*4*pi*R3/lambda) ; 

s3o  =  (shiftM(s3u’  , -2*shift_incr) )  ’  ;  ”/0coregistration. 
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"/compute  fractional  Fourier  transforms 

alpha  =  acot (k*vrel~2*N/ (pi*Rb*f s~2) ) ;  ‘/calculate  theoretical  alpha 
a  =  alpha/ (pi/2) ; 
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dpcaMag234(n)  =  max(abs(DPCA234)) ; 
DPCA234_clut  =  dpca23_clut+dpca34_clut ; 
dpcaMag234_clut (n)  =  max(abs(DPCA234_clut)) ; 
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Annex  E:  MATLAB  Code  ‘airMultiChan  toggle. m 
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°/0Target  signal  length  Ts: 

Theta  =  3;  °/„3  degree  for  3dB  beamwidth 


theta  =  Theta*pi/180 ;  "/convert  to  radians 

az  =  int32 (Rb*tan(theta/2) *5 . 14) ;  "/convert  beam  angle  to  sample  number 
az  =  double  (az);  '/[samples] 

Ts  =  2*az*tpri;  °/„[s]  target  signal  length 
T  =  3.005*Ts;  "/total  signal  length 
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sf  =  sum(sl ( : ,ml) ,2) ; 

scl  =  sum(sl ( : ,m2) ,2) ;  sclr  =  shiftM(scl,-((8*D/16)/va)/tpri) ; 

sc2  =  sum(s2( : ,ml) ,2) ;  sc2r  =  shiftM(sc2,-((8*D/16)/va)/tpri  +  1/2);  %sclr  =  sc2r; 
sa  =  sum(s2 ( : ,m2) , 2) ;  sar  =  shiftM(sa, - ( (16*D/16) /va) /tpri  +  1/2); 
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sa_clut  =  rectpuls (t2 (tb+delta_t) ,Ts) . *sum(s2_clut ( : ,m2) ,2) ;  sa_clutr  =  shiftM(sa_clut , - ( (16*D/16) /va) /tpri+1/2) ; 
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theta  =  Theta*pi/180 ;  "/convert  to  radians 

az  =  int32 (Rb*tan(theta/2) *5 . 14) ;  "/convert  beam  angle  to  sample  number 
az  =  double  (az);  '/[samples] 

Ts  =  2*az*tpri;  °/„[s]  target  signal  length 
T  =  3.005*Ts;  "/total  signal  length 
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sc  =  sum(s2( :  ,m2)  ,2) ;  scr  =  shiftM(sc,-((4*D/16)/va)/tpri+l/2) ;  '/need  to  modify  this 
sa  =  sum(sl ( : ,m3) , 2) ;  sar  =  shiftM(sa, - ( (8*D/16) /va) /tpri) ; 
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"/calculate  theoretical  fractional  angle 
N  =  length (tl); 

frac_a  =  (2/pi) *acot (beta*N*vrel~2/ (pi*Rb*f s~2) ) ; 
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Annex  G:  MATLAB  Code  ‘airMultiChan  fourChan.m’ 
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°/0Target  signal  length  Ts: 

Theta  =  3;  °/„3  degree  for  3dB  beamwidth 


>■> 

p 

o 

-p 

0 

•H 

A 

0 

bO 

p 

0 

p 

!>> 

P  CD 
O  > 
-P  -H 
CO  CD 
•H  CJ 
rP  CD 
P 
CD 
bO 
P 


0 

P 

P 

CD 

CN 

P 

P 

< 

.  ~ 

P 

•H 

rP 

p 

B 

P 

0 

+ 

T— 1 

p 

P 

CD 

0 

CN 

fP 

i — 1 

P 

< 

Ph 

P 

+ 

P 

o-£ 

P 

Co 

0 

H — 1 

CO 

.  ~ 

B 

P 

l 

Ptf 

o 

T— 1 

' - - 

p 

CN 

* 

< 

0 

CD 

■P 

+ 

P 

i — 1 

0 

bO 

+ 

CN 

P 

P 

< 

* 

cd 

CN 

•r— ) 

< 

1 

P 

' — ^ 

cd 

P 

CO 

Oh 

0 

0 

0 

H — 1 

X 

P 

0 

B 

\ 

0 

0 

l 

Q 

* 

P 

i — 1 

T — 1 

* 

CO 

P 

rP 

P 

0 

P 

p 

0 

P 

Ph 

0 

0 

0 

0 

> 

bO 

0 

+ 

B 

H 

•H 

P 

P 

X3 

i — 1 

P 

l 

x* 

o 

0 

X* 

P 

O 

CN 

H — 1 

0 

o 

i— 1 

O 

Ph 

P 

< 

B 

P 

P 

0 

•  ~ 

1 

i— 1 

rP 

P 

•H 

.  *N 

/~N 

l 

0 

o 

.  *N 

0 

P 

. « 

0 

rP 

t— 1 

/— N 

S  N 

P 

p 

P 

bO  CN 

> 

1 - 1 

P 

P 

r— 1 

bO 

P 

\ 

0 

0 

B 

H - 1 

0 

0 

0 

p 

H — 1 

•H 

0 

o-£ 

bO 

' - - 

B 

B 

X* 

p 

0 

I- 1 

•H 

P 

rP 

l 

l 

CD 

LO 

P 

.  «N 

0 

P 

H - 1 

T — 1 

H - 1 

+ 

> 

* 

r~ -i 

P 

I- 1 

Oh  CN 

P 

bO  -p 

P 

P 

P 

S  N 

CO 

0 

0 

p 

N 

P 

N — S 

* 

* 

P 

o 

CN 

0 

bO 

P 

1 

H 

i — 1 

0 

,P 

.  *N 

.  ~ 

0 

0 

P 

o 

\ 

i — 1 

P 

b0  H 

0 

1 — 1 

P 

P 

p 

> 

> 

' — ^ 

0 

Ph 

0 

•H 

' — ' 

•H 

P 

bO 

0 

0 

_ _ • 

' — ' 

P 

c^ 

P 

0 

P 

bO 

H — 1 

P 

B 

B 

1 

. « 

0 

•H 

Ph 

•H 

1 _ 1 

0 

l 

1 

1 

1 

o 

,0 

0 

f— 1 

i— 1 

p 

P 

0 

i — 1 

t — i 

T — 1 

P 

00 

P 

i _ i 

0 

0 

Ph 

i _ i 

P 

P 

P 

P 

0 

T — 1 

i _ i 

P 

P 

CN 

P 

H - 1 

* 

* 

0 

0 

B 

\ 

P 

O 

\ 

•H 

P 

B 

o 

o 

B 

B 

l 

•H 

0 

.  « 

P 

CN 

/— N 

B 

' - - 

n — y 

X 

1 

l 

H - 1 

Oh  P 

.  « 

o-£ 

\ 

•H 

0 

P 

P 

> 

> 

T— 1 

H — 1 

P 

* 

* 

N 

•H 

P 

P 

0 

0 

X 

X 

' - - 

0 

P 

0 

P 

.  « 

•H 

Ph 

0 

.  ~ 

B 

B 

+ 

+ 

_ _ ' 

' — ^ 

0 

P 

pZ 

' — ^ 

Ph 

0 

P 

P 

P 

LO 

Ph 

Ph 

i — 1 

CD 

0 

P 

H 

Ph 

* 

P 

T— 1 

0 

0 

o 

o 

_ _ ' 

' — ^ 

P 

rP 

CN 

i — 1 

* 

* 

P 

CN 

P 

P 

X 

!>b  P 

p 

Oh 

H 

CO 

P 

N 

LO 

1 

1 

0 

CN 

p 

p 

P 

P 

P 

0 

o 

H 

H 

P 

II 

II 

II 

II 

cr1 

CJ* 

O 

II 

P 

o 

* 

o 

' — ^ 

n — y 

0 

LO 

0 

0 

0 

•H 

xJ 

CN 

1 

1 

P 

T — 1 

P 

P 

p 

P 

P 

cd 

CO 

0 

1 

0 

0 

0 

0 

ll 

II 

P 

II 

II 

II 

II 

II 

P 

B 

B 

B 

B 

II 

CD 

II 

0 

II 

1 

1 

1 

l 

T — 1 

t — 1 

,P 

N 

N 

0 

H - 1 

CN 

bO 

■H 

H — 1 

H - 1 

T — 1 

P 

P 

T— 1 

P 

0 

0 

E-< 

H 

P 

P 

B 

P 

s 

X 

t>>  P5 

P 

0 

52 


>> 

u 

o 

p 

CO 

•H 

d 

a) 

bO 

d 


d 

P 

a) 

o 

> 

p 

•H 

0 

0 

•H 

(J 

,0 

0 

P 

0 

bO 

P 

.  ~ 

0 

s — % 

P 

CN 

P 

< 

.  ~ 

•H 

rP 

s — s 

B 

0 

+ 

CN 

p 

P 

0 

CN 

CP 

P 

< 

P 

+ 

o-£ 

P 

0 

CN 

. « 

B 

P 

/•— N 

l 

CP 

CN 

' — ^ 

CN 

* 

< 

0 

,P 

+ 

P 

0 

+ 

CN 

rP 

< 

* 

CN 

< 

’I-) 

1 

' - - 

P 

CO 

Ph 

0 

T — 1 

X 

B 

\ 

0 

l 

Q 

* 

CN 

* 

P 

f — N 

0 

0 

+ 

B 

H 

l 

CN 

CN 

< 

B 

p 

1 

AC 

S'  N 

1 

0 

i~ n 

/^~N 

P 

ss. 

(— 1 

P 

p 

r— 1 

B 

T — 1 

0 

0 

0 

' — ^ 

B 

B 

xJ 

rP 

/~N 

l 

l 

P 

CN 

CN 

CN 

+ 

bO  -p 

P 

P 

P 

_ _ ' 

* 

* 

rO 

0 

0 

0 

P 

1 - 1 

P 

P 

p  > 

> 

' — ^ 

P 

h0 

0 

cb  ^ 

P 

0 

H — 1 

P 

a 

S 

1 

0 

0 

.  *N 

.  ~ 

. « 

1 _ 1 

0 

1 

1  1 

1 

e 

p 

1 - 1 

CN 

CN 

p 

P 

p 

CN 

CN 

CN 

1 _ 1 

P 

P  P 

P 

0 

h0 

p 

CN 

* 

* 

*  O 

0 

a 

•rH 

d 

/^N 

P 

6 

O 

o  g 

a 

1 

0 

0 

H - 1 

CN 

CO 

X 

1 

1 

CN 

0 

Oh 

s 

B 

B 

P 

P 

> 

>  CM 

CN 

P 

0 

0 

0 

0 

X 

X 

1 

H 

H 

+ 

+  ^ 

' — ^ 

0 

0 

o 

' — ^ 

Oh 

Oh 

i — 1 

P 

p 

H - 1 

CN 

H — 1 

0 

0 

o 

o  ^ 

' — ^ 

P 

P 

.  ~ 

0 

0 

0 

0 

d 

d 

X 

p 

p 

Oh 

P 

.  ~ 

CO 

d 

d 

P 

d 

.  «N 

.  «N 

CN 

H — 1 

0 

B 

s 

B 

II 

II 

II 

II  cr1 

cr1 

o 

0 

00 

H — 1 

p 

P 

P 

P 

0 

0 

0 

Oh 

CO 

0 

0 

0 

0 

P 

P 

p 

P 

d 

0 

H - 1 

LO 

CT) 

T — 1 

d 

0 

0 

0 

0  II 

ll 

1 

0 

II 

II 

II 

s 

6 

6 

s 

II 

o 

II 

II 

II 

II 

P 

1 

1 

1 

1  CN 

CN 

P 

0 

H - 1 

CN 

CO 

CN 

CN 

CN 

CN  P 

d 

CN 

0 

H - 1 

CN 

CO 

bO  d 

d 

P 

u 

X 

0 

s 

B 

B 

B 

u 

u 

o 

DRDC  Ottawa  TM  2006-209 


p 

O 

pZ 
+ 
t — i 

P 

u 

PZ 

* 

P 

-p 

CD 

P 

* 


CN 

P 

U 

CP 

+ 

CN 

-P 

U 

' — ^ 
* 
P 
-P 
CD 
P 
* 


. « 

Ph 

. ~ 

Ph 

^~N 

CN 

X 

/~\ 

X 

CN 

\ 

CN 

CD 

CN 

CD 

\ 

T — 1 

< 

* 

< 

* 

.  « 

. « 

. « 

. « 

T— 1 

,P 

A 

/— N 

^~N 

/  \ 

+ 

/*  N 

s  \ 

CN 

CN 

CN 

CN 

+ 

.  »N 

+ 

0 

+ 

CO 

* 

•* 

* 

•* 

/— n  *H 

O  P 
K5’"!  I 

P 


o  -p 
>>  i 
cd 


-rH  CN  00  ^ 

a  a  a  a 


•H  -H  p 

P  P  Oi 

Ph  Ph  P 
P  P  \ 

\  \  /— N 

cd 


+ 

p 

+ 

p 

P 

p 

P 

p 

p 

p 

> 

. « 

rH 

.  ~ 

rH 

P 

p 

P 

p 

> 

> 

\ 

f  N 

CN 

0 

CN 

0 

rH 

rH 

i — 1 

rH 

\ 

\ 

CN 

< 

Xl 

CN 

< 

X3 

O 

U 

o 

o 

CD 

< 

+ 

< 

+ 

1 

1 

1 

1 

CD 

CD 

T — 1 

.  *N 

nP 

p 

rP 

/~N 

P 

T— 1 

CN 

H — 1 

CN 

H — 1 

T — 1 

\ 

^~N 

p 

/'-'N 

P 

0 

0 

0 

0 

\ 

\ 

Q 

+ 

CD 

rH 

+ 

CD 

rH 

- - ' 

_ _ ' 

Q 

Q 

* 

CN 

H — 1 

O 

T— 1 

O 

g 

g 

s 

s 

* 

* 

CN 

< 

CN 

\ 

1 

CN 

\ 

1 

p 

p 

p 

p 

00 

H - 1 

0 

< 

Q 

P 

< 

Q 

P 

0 

0 

0 

0 

0 

' — ^ 

P 

O 

* 

P 

O 

* 

P 

i — 1 

* 

* 

* 

* 

' — ^ 

w 

V—/ 

* 

>>  P 

_ _ ' 

>>  P 

' — s 

P 

1 

1 

1 

P 

P 

1 

P 

1 

P 

/^~N 

/'—N 

s 

•V 

•v 

CP 

+ 

a 

i 

P 

+ 

a 

i 

P 

b0 

0 

0 

0 

0 

CN 

CO 

0 

* 

P 

P 

•H 

H 

H 

H 

H 

,P 

,p 

,P 

rH 

•H 

CN 

H — 1 

a 

CN 

CN 

a 

0 

o 

o 

u 

bO 

Ph 

< 

a 

i 

< 

a 

i 

S  N 

/"  N 

/^~N 

N 

' — ^ 

Vw/ 

w 

P 

' — ^ 

H - 1 

CN 

P 

P 

P 

P 

P 

S 

s 

S 

P 

\ 

f  N 

i 

P 

/—■ \ 

i 

P 

0 

1 

1 

1 

1 

p 

p 

p 

CN 

/■  N 

f  \ 

' — s 

P 

P 

P 

P 

P 

HH 

4H 

MH 

rH 

< 

P 

p 

0 

p 

p 

0 

P 

P 

P 

P 

P 

•H 

•H 

•H 

P 

i — 1 

0 

P 

p 

i — i 

p 

p 

i — i 

P 

I— 1 

I— 1 

rH 

rH 

0 

rP 

x\ 

rP 

P 

0 

0 

a 

i 

a 

i 

P 

a 

i 

a 

i 

P 

r— 1 

0 

0 

0 

0 

i — 1 

0 

0 

0 

O 

P 

P 

Ph 

Ph 

o 

xJ 

Td 

X* 

P 

•H 

> 

,P 

T— 1 

■H 

P 

CN 

CN 

P 

+ 

+ 

+ 

+ 

P 

II 

II 

II 

P 

* 

Ph 

.  *N 

P 

P 

O 

P 

p 

O 

p 

P 

P 

P 

P 

bO 

O 

S 

P 

* 

* 

0 

* 

* 

0 

o 

P 

P 

P 

P 

•H 

P 

P 

P 

P 

* 

1— 1 

. « 

CO 

> 

P 

P 

P 

P 

p 

P 

MH 

_ _ ' 

_ _ ' 

' - - 

0 

CN 

CO 

P 

P 

H 

/  N 

1 — 1 

* 

> 

> 

> 

> 

1 

1 

1 

1 

rP 

,P 

nP 

MH 

P 

< 

P 

cd 

- - ' 

' — / 

+ 

w 

_ _ ' 

+ 

0 

P 

U 

u 

u 

0 

> 

. « 

P 

CN 

0 

H - 1 

CN 

H - 1 

CN 

0 

rH 

rP 

i — 1 

1 

b0\ 

1 

1 

P 

1 

1 

P 

P 

P 

P 

P 

P 

P 

.  *N 

.  -N 

. « 

. « 

P 

' - - 

P 

o 

CN 

•H 

H 

. « 

P 

P 

P 

_ _ ' 

_ _ ' 

' - - 

P 

P 

P 

p 

p 

O 

P 

u 

X 

CO 

P 

o 

u 

rH 

O 

U 

rH 

P 

0 

0 

0 

0 

P 

P 

P 

p 

p 

•H 

o 

•H 

> 

/■— \ 

•H 

> 

X 

X 

o 

X 

X 

U 

P 

rH 

rH 

i — 1 

i — 1 

rH 

I- 1 

rH 

rH 

rH 

P 

o 

P 

_ _ ' 

■n^ 

p 

P 

\ 

_ _ ' 

' — .✓ 

1 

' — ■/ 

_ _ ' 

1 

0 

P 

P 

P 

P 

O 

u 

U 

u 

o 

0 

p 

0 

B 

CD 

Ph 

o 

_ _ ' 

' — ✓ 

H - 1 

' — s 

- - ' 

CN 

Ph 

Ph 

Ph 

Ph 

Ph 

1 

1 

1 

1 

P 

* 

P 

+ 

P 

P 

X 

p 

p 

CO 

p 

p 

0 

P 

P 

P 

P 

P 

X5 

T— 1 

CN 

CO 

O 

.  -N 

✓-N 

O 

p 

P 

* 

/~N 

p 

p 

p 

p 

1 

U 

O 

O 

o 

P 

rP 

rP 

rP 

0 

s~\ 

•H 

0 

' - S 

P 

o 

P 

II 

cr1 

cP 

II 

cr1 

cr1 

II 

P 

0 

0 

0 

0 

P 

U 

u 

u 

u 

rP 

H — 1 

Ph 

rP 

> 

CN 

r— 1 

LO 

X 

.  ~ 

.  «N 

0 

0 

0 

CO 

P 

P 

P 

P 

p 

P 

P 

\ 

P 

\ 

CO 

o 

_ _ ' 

O 

o 

P 

p 

p 

T — 1 

P 

0 

P 

+ 

+ 

+ 

+ 

' - - 

CN 

CN 

w 

CN 

nP 

X 

p 

ll 

II 

p 

1 

ll 

ll 

P 

II 

II 

II 

II 

0 

0 

rP 

' — ^ 

0 

< 

a 

CD 

\ 

P 

II 

II 

I— 1 

rH 

P 

rH 

0 

bO  ■h 

CN 

CO 

P 

P 

P 

rH 

p 

P 

H 

bO 

II 

u 

T— 1 

H — 1 

o 

CN 

CN 

O 

P 

P 

P 

P 

p 

P 

rP 

rP 

nP 

A 

P 

bO 

II 

P 

0 

CO 

P 

1 

P 

P 

p 

1 

P 

p 

1 

II 

p 

P 

1 

P 

P 

P 

P 

p 

P 

o 

u 

o 

u 

rH 

P 

1 - 1 

P 

P 

_ _ ' 

0 

P 

p 

u 

p 

o 

o 

H — 1 

o 

o 

CN 

P 

rH 

rH 

rH 

rH 

P 

P 

0 

P 

P 

> 

ii 

CD 

rH 

rH 

rH 

X 

p 

CP 

0 

P 

CP 

CP 

0 

0 

O 

U 

O 

o 

II 

II 

II 

II 

O 

i — 1 

1 

(J 

P 

II 

U 

u 

P 

1 

1 

1 

1 

xJ 

rH 

U 

1 - 1 

II 

CD 

II 

1 

1 

p 

X* 

0 

T — 1 

CN 

CO 

Td 

■H — 1 

CN 

CO 

P 

II 

P 

P 

rP 

h0 

P 

H — 1 

CN 

o 

P 

bO  .d 

rP 

,P 

P 

,P 

rP 

rP 

,p 

O 

P 

o 

0 

u 

X 

P 

CO 

CO 

MH 

0 

o 

U 

o 

u 

u 

u 

u 

u 

S 

P 

> 

DRDC  Ottawa  TM  2006-209 


53 


atiPhase  =  (beta* (D/2) *gamma*vyO/va  +  beta* (D/2) *ve*delta_t/Rb) *180/pi  °/0toggle  mode  baseline  is  D  not  D/2 
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atiPhase_max  =  angle (ati2 (J, 1) ) *180/pi ; 

title(  [’Measured  ATI  Phase  2  &  3  =  num2str (atiPhase_max)] ) ; 
xlabel( [’Theoretical  ATI  Phase  =  ’,  num2str (atiPhase)] ) ; 
ylabel([’ATI  Phase  Error  =  ’,  num2str (abs (atiPhase-atiPhase_max) )] ) ; 
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